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Two-dimensional eigenvalue analysis is used on a massive scale to study spatial 
instabilities of compressible shear flows with two inhomogeneous directions. The main focus 
of the study is crossflow dominated swept-wing boundary layers although the methodology 
can also be applied to study other type of flows, such as the attachment-line flow. Certain 
unique aspects of formulating a spatial, two-dimensional eigenvalue problem for the 
secondary instability of finite amplitude crossflow vortices are discussed, namely, fixing the 
spatial growth direction unambiguously through a non-orthogonal formulation of the 
linearized disturbance equations. A primary test case used for parameter study corresponds 
to the low-speed, NLF-0415(b) airfoil configuration as tested in the ASU Unsteady Wind 
Tunnel, wherein a span wise periodic array of roughness elements was placed near the 
leading edge in order to excite stationary crossflow modes with a specified fundamental 
wavelength. The two classes of flow conditions selected for this analysis include those for 
which the roughness array spacing corresponds to either the naturally dominant crossflow 
wavelength, or a subcritical wavelength that serves to reduce the growth of the naturally 
excited dominant crossflow modes. Numerical predictions are compared with the measured 
database, both as indirect validation for the spatial instability analysis and to provide a basis 
for comparison with a higher Reynolds number, supersonic swept-wing configuration. 
Application of the eigenvalue analysis to the supersonic configuration reveals that a broad 
spectrum of stationary crossflow modes can sustain sufficiently strong secondary instabilities 
as to potentially cause transition over this configuration. Implications of this finding for 
transition control in swept wing boundary layers are examined. 


Nomenclature 


B,C 


N y 

U s 


= initial amplitude of stationary crossflow vortex 
= coefficient matrices of the linearized Navier-Stokes equations 
= N factor 
= time 

= number of points used in the wall normal direction in the numerical discretization 
= number of points used in the direction parallel to the wing leading edge in the numerical discretization 
= velocity in the direction aligned with local inviscid streamline 


V, 


subscripts 
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y 

z 

X 
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coefficient matrices of the linearized Navier-Stokes equations 

wing surface coordinate in the direction perpendicular to the leading edge 
coordinate normal to the wing surface 

wing surface coordinate in the direction parallel to the leading edge 
wing surface coordinate along a crossflow vortex line 
wing surface coordinate perpendicular to a crossflow vortex line 
complex wave number 
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= real part of OC 
= imaginary part of CC 

= angle between a crossflow vortex line and a surface line perpendicular to the leading edge. 
= angular frequency 

= wing surface coordinate along a crossflow vortex line in a non-orthogonal system 
= wing surface coordinate parallel to the leading edge in a non-orthogonal system 


I. Introduction 

Even in a low-disturbance environment that is characteristic of external aeronautical flows, laminar-turbulent 
transition in a 3D boundary layer can occur through many routes. Potential paths for transition include the 
destabilization of attachment line boundary-layer, the development and breakdown of streamwise instabilities (i.e., 
Tollmien-Schlichting waves and the first and second modes of instability depending on flow Mach number), and the 
growth and breakdown of crossflow vortex instabilities (if the boundary layer is three-dimensional) or the 
centrifugal or Goertler mode instability (if the surface geometry includes regions of concave streamwise curvature) 
[1, 2, 3]. For flow configurations that are dominated by one of the streamwise instabilities, efficient and usually 
reliable predictions of transition onset can be achieved via linear growth (N-factor) correlations based on classical 
stability theory or linear PSE (Parabolized Stability Equations) [4, 5]. However, there is increasing experimental 
evidence that receptivity and nonlinear processes become relatively more significant in the case of transition due to 
crossflow or Goertler vortex instabilities [6-9]. A more holistic prediction approach based on the nonlinear 
evolution of these vortices and the ensuing high-frequency secondary instabilities then becomes desirable [10-12]. 

An earlier application of a higher-fidelity transition prediction approach [10] had yielded promising results in the 
context of a low-speed crossflow experiment involving a 45-deg swept NLF-0415(b) airfoil [6]. Specifically, an re- 
factor correlation based on the linear amplification of secondary instabilities of stationary crossflow vortices of 
known amplitude provided a more robust correlation with the measured location for transition onset than an absolute 
amplitude correlation based on the primary instability alone. The first known application of secondary N-factors in 
the context of semi-empirical transition prediction was presented by El-Hady in the context of streamwise 
instabilities in 2D supersonic boundary layers (see, for instance, Ref. [13]). The preliminary results based on this 
higher fidelity approach were also consistent with the observed delay of transition in the presence of a subcritical 
stationary crossflow mode that was seeded via artificial, spanwise periodic roughness near the leading edge. 

The technique of crossflow-transition control via spanwise periodic roughness is also relevant to supersonic swept- 
wing flows, but has not yet been demonstrated beyond initial laboratory experiments in the Mach 2.4 Supersonic 
Wind Tunnel at Arizona State University (ASU) [14]. A preliminary application of the higher fidelity transition 
prediction approach to the latter experiments was presented in Ref. [11, 12]; however, it did not include an explicit 
account of the secondary instability phase during transition. A major goal behind the present work is to extend the 
findings in Ref. [11] by including the growth of high-frequency secondary instabilities of stationary crossflow 
modes and, more generally, to further develop the holistic prediction approach in the context of compressible 
boundary layers. To that end, the temporal instability analysis of Refs. [10, 15] is extended to spatial instability 
modes, allowing systematic mode tracking without any additional simplifying assumptions as required during the 
temporal analysis. It should be noted that such spatial predictions can be directly compared with numerical 
simulations and, perhaps, used to provide suitable initial conditions for a simulation that is focused on the (strongly 
nonlinear) laminar breakdown phase. 

Mathematically, prediction of stationary crossflow transition via the secondary N-factor approach differs from the 
classical N-factor analysis based on (small-amplitude) primary instability modes, in that the basic state for 
secondary instability analysis (which has been modified by the presence of finite amplitude stationary vortex 
instabilities) is strongly inhomogeneous in two spatial directions (wall normal and spanwise) instead of just a single 
direction (wall normal) as during the classical analysis using the unperturbed boundary layer flow as the basic state. 
Consequently, the amplification characteristics and mode shapes of the high-frequency secondary instability modes 
are governed by an eigenvalue problem based on a two-dimensional partial differential equation (PDE) [3, 10, 16, 
17], rather than the simpler, classical eigenvalue problem that is derived from a system of ordinary differential 
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equations (ODEs). Whereas the latter set of ODEs involves the wall-normal coordinate as the only independent 
variable, the 2D eigenvalue problem for the secondary instability of finite amplitude crossflow modes involves an 
additional, surface-tangential coordinate that cuts across the cross section of the crossflow vortex. In cases where 
the vortex pattern involves a single dominant wavelength, the secondary perturbations may be treated as being 
periodic along this latter coordinate. Attachment line flows are not considered in the present paper, but it should be 
noted that linear (primary) instabilities of an attachment line boundary layer can also be described using a 2D 
eigenvalue problem [18-20], except that both the basic state and the unstable perturbations now have an algebraic, 
rather than periodic, behavior along the second spatial coordinate. 

A number of hydrodynamic instability analyses based on partial differential eigenvalue problems have been reported 
in the literature (see [3, 10, 16-21] and the references therein), mostly in the context of the temporal instability (such 
that the perturbations are spatially periodic and amplify in time) and/or incompressible flows, with only a few 
analyzed compressible flows (e.g. [16, 19]). The present paper is focused on the development and application of 
spatial instability analyses for subsonic and supersonic swept-wing boundary layers. Spatially growing secondary 
instabilities of crossflow vortices over an infinite-span airfoil have been considered. The extension from temporal to 
spatial predictions is not quite straightforward for secondary instabilities of crossflow vortices; the additional 
complexity arises because of the differences in the orientation of the vortex axis (typically within a few degrees of 
the inviscid streamline) along which the basic state evolves on a slow spatial scale and the direction in which the 
secondary instability modes amplify (viz., along the chordwise direction). A brief summary of the unique aspects 
of spatial secondary instability analysis for infinite span swept-airfoil flows is presented in section II below. Section 
III presents numerical results pertaining to the low-speed, NLF-0415(b) configuration and compare the predictions 
with measured data, both as a qualitative validation for the spatial instability predictions and to provide a reference 
to assess similar results for a supersonic configuration modeling an experiment in the ASU Mach 2.4 Supersonic 
Wind Tunnel (Section IV) [14]. Concluding remarks are presented in section V. 


II. Non-orthogonal Formulation of Secondary Instability Equations and Solution Precedure. 

The equations governing the unstable perturbations are obtained by linearizing the Navier-Stokes equations about a 
specified basic state, e.g., a finite amplitude crossflow vortex developing in a swept-airfoil boundary layer. For a 
basic state that is slowly varying along one spatial coordinate (viz., the vortex axis in the case of crossflow vortices), 
using a wave ansatz in that direction for the disturbance quantities (in the spirit of WKB or multiple scale analysis) 
reduces the spatial dimension of the problem by one, resulting in a set of two-dimensional, linear partial differential 
equations at the leading order. In the case of temporal secondary instability analysis for stationary crossflow 
vortices in a swept airfoil boundary layer [3, 10, 16], the coordinate system is usually rotated so that one of the 
coordinate axes along the surface is aligned with the local vortex direction ( X ) and the other surface coordinate 
(Z) is orthogonal to X (Fig. 2.1). This leads to a set of linear partial differential equations in terms of the surface 
normal coordinate Y and the across-the-vortex coordinate Z . Together with appropriate boundary conditions at the 
wall and free-stream boundaries, and periodic boundary conditions in Z , one is led to a planar (i.e., 2D) eigenvalue 
problem which may be solved using numerical techniques for large linear systems. 

The above choice of surface coordinates becomes inappropriate in the context of spatial instability analysis. 
Specifically, the spatial exponential growth obtained from the 2-D eigenvalue problem in the coordinate system 
(X,Z) is in the direction of the vortex axis and, therefore, has a component in the direction parallel to the leading 
edge, which is not the desired solution. In addition, the base flow is not strictly periodic in Z, so that approximations 
have to be made based on the fact that the deviation from periodicity is typically small. A better approach is to seek 
a solution that is consistent with the characteristics of the underlying basic state. Specifically, for a basic state 
involving crossflow vortices of fixed spanwise wavelength, one needs to seek secondary instability modes that are 
periodic in the direction parallel to the leading edge ( £ or Z - direction in Figure 2.1) and grow only in the direction 
perpendicular to the leading edge ( X -direction in Figure 2.1). 

For infinite-span configurations, the unperturbed boundary layer is invariant in the direction z parallel to the 
leading edge (Fig. 2.1), which allows both the primary (i.e., the stationary crossflow vortex) and the secondary 
disturbance fields to be periodic in that direction. However, because the primary disturbance field exhibits “rapid” 
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variation in the chordwise coordinate X , the orthogonal coordinate system (x, Z.) cannot be used to reduce the 
spatial secondary instability analysis to a 2D eigenvalue problem. However, the latter goal can be achieved by using 
a non-orthogonal, vortex aligned coordinate system E, and E, (where E, = X and E, = z ) as described below. It is 
worth noting that the non-orthogonal formulation precisely mimics the DNS computations in refs. [22-23], wherein 
the vortex aligned non-orthogonal coordinate system was employed towards efficient simulation of disturbance 
evolution in swept airfoil boundary layers without any approximations. For general 3D configurations, the 
implementaton of parabolized stability equations (PSE) in the LASTRAC code [24] also allows the option of a non- 
orthogonal coordinate system. 

The linearized compressible Navier-Stokes equations in the non-orthogonal system can be written in the following 
form 


dt df 
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where (f) is a vector whose elements are the perturbation pressure, velocities and temperature, the coefficients of the 
derivatives are 5x5 matrices that are slow varying functions of E, , functions of y and periodic functions of E, . 

In the surface coordinate system ( < ^ , C, ) , the base flow is slowly varying in the E, -direction and periodic in the E, - 
direction. The perturbation variables associated with the secondary instability can then be expressed in the form 

0(£ y,C) = v(€>y>Chx p(- icot + i { ad%) (2) 

where y is the coordinate normal to the wing surface, <p{^ , V, ^ ) is the 2-D complex eigenfunction at station E, 
that is periodic in E, ; CO is the given frequency; and (X is the complex wave number (i.e., the spatial eigenvalue). 

It should be emphasized that, in the non-orthogonal coordinate system, the use of CXd^ in the exponential part does 
not constitute allowing the disturbance to grow in the spanwise (z) direction. This point will become clearer in the 
course of later discussion, which presents the expression for the evolution of perturbation quantities in the more 
intuitive orthogonal coordinate system. 

On substituting Eq. (2) into Eq. (1) and neglecting the dependencies on E, of both the coefficients of the equations 
and the mode shape (p , one obtains 
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Equation (3) is solved as a 2-D spatial eigenvalue problem with OC as the eigenvalue for a given CO . Periodic 
boundary condition is used in the E, direction and homogeneous Dirichlet conditions are used for the velocity and 
temperature perturbations at the wall, along with suitable boundary conditions along the free-stream boundary. 
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The form of the solution in Eq. (2) is not intuitive and, for better clarity, one can rewrite it in terms of the more 
intuitive orthogonal coordinates. To this end, it is easy to verify that the surface coordinate systems ( X , Z.) and 
(g, are related through the following expressions 


d% 


dx 

cos# 


(4) 


= dz- dx tan# 


(5) 


where # denotes the angle between the g and the X axes. Therefore, the right hand side of Eq. (2) can be written 

f r dx 
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y J cos# J cos # J 


where CC r and OC i are, respectively, the 


real and imaginary parts of Ci . Observe that, at any fixed X , (p is a periodic function of Z ; the solution oscillates in 
X along with an exponential growth or decay in its amplitude; and, furthermore, there is no exponential growth in z. 
The N-factor is then defined as 

AT r adx 

N =-j- L -z (6) 

J cos# 

where the lower limit of integration corresponds to the onset of the instability. 


The formulation is also valid when the base flow and, hence, the coefficient matrices in (1) are not periodic 
in £ direction. In this case, the eigenfunction <p{g , V, d, ) is also not periodic in d, and an alterative set of boundary 
conditions needs to be applied in that direction. 


Numerically, Equation (3) is discretized in the y-direction using high order finite difference and, when the 
coefficient matrices are periodic, discrete Fourier spectral method is used in the C, -direction, otherwise high order 

finite difference is used. The resulting algebraic eigenvalue problem has a dimension of 5 X Nr X N , however, its 

banded structure can be utilized for storage and computational efficiency. The algorithms used for solving the 
algebraic eigenvalue prolem are the same as those used in previous analyses of 2D eigenvalue problems (e.g. [10]). 


A typical N-factor computation procedure involves several steps: 


(0) Computation of basic state for secondary instability: The unperturbed boundary-layer flow is computed by 
solving compressible boundary-layer equations for a specified inviscid pressure distribution that is derived 
from Euler calculations. The evolution of stationary crossflow vortices in this flow is computed using 
nonlinear Parabolized Stability Equations (PSE) as implemented in the Langley Stability and Transition 
Analysis Code (LASTRAC) [24], The wall-normal resolution used for these computations varied between 
121 to 251 points, and between 64 to 96 Fourier modes were typically used in the span wise direction. For 
secondary instability calculations, the basic state computed with nonlinear PSE was interpolated onto a grid 
that is more suitable for secondary instability modes. A typical resolution used for the eigenvalue 
computations presented in this paper corresponds to 121 points in the wall-normal direction and 32 points 
in the periodic ( d, ) direction, although additional resolution was found to be necessary in some of the 
cases. The computational grid was tailored to the anticipated/computed mode shapes and spot checks were 
employed to verify that the shift in eigenvalues was insignificant when the number of points was further 
increased and/or when the clustering associated with the wall-normal grid was modified. 

(1) A thorough eigenvalue search is performed at one or more chordwise locations for a range of 
wavenumbers. The Arnoldi method is used for global computations of temporal eigenvalues which, in turn, 
are used as starting values for an iterative spatial computation. This establishes the number of unstable 
modes and the corresponding range of frequencies. 

(2) A suitable subset of eigenvalues for each mode is chosen to span the relevant range of frequencies. This 
forms the starting set of eigenvalues to be used for N-factor calculations for that mode. 
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(3) With each selected eigenvalue from step (2) as a starting guess, the eigenvalue computation is marched 
both upstream and downstream in the chordwise direction to cover the appropriate range of locations. At 
each step during the marching process, the starting guess is updated via linear extrapolation from the 
previously computed eigenvalues at adjacent locations. Having determined the eigenvalues over the region 
of interest, the N-factor evolution for each frequency and mode type is calculated using Eq. (6). 

III. Secondary Instability of Crossflow Vortices: Low Mach Number NLF-0415b Case 

This section examines the secondary instability of crossflow vortices for the NLF-0415b configuration that was used 
for detailed measurements of crossflow induced transition in the Arizona State University (ASU) Unsteady Wind 
Tunnel [25]. During those experiments, a span wise periodic array of roughness elements was placed near the airfoil 
leading edge to introduce stationary crossflow vortices with the desired wavelength and varying initial amplitude. 
Two sets of experimental data points are modeled herein, namely those corresponding to the vortex wavelengths of 
12 mm (the naturally dominant stationary crossflow mode) and 8 mm (subdominant mode that was found to weaken 
the growth of the naturally excited 12 mm mode and, hence, to delay the onset of transition on the airfoil), 
respectively, at the fixed chord Reynolds number of 2.4 million. Similar to [10], initial amplitudes of the stationary 
vortex are chosen so that the computational amplitudes at a selected upstream location are closely matched with the 
experimental results for three roughness heights of 6, 18 and 48 microns, respectively, for the 12 mm case. Although 
the roughness induced receptivity is not directly modeled in this paper, for convenience of notation, these cases will 
be referred to on the basis of the corresponding experimental configuration. Accordingly, the above three cases are 
denoted as H06S12, H18S12 and H48S12, respectively. Similarly, for the 8 mm case, two roughness heights of 6 
and 48 microns are chosen to match the computational results. These are referred to as H06S08 and H48S08, 
respectively. The stabilizing effect of a higher amplitude 8 mm mode on the 12 mm stationary mode has already 
been established in both experiments [27] and computations [10]. Thus, the goal behind the X z = 8 mm 
computations presented in this paper is to examine the growth potential for secondary instabilities at that wavelength 
and their likely effect on the transition onset location. 

In Fig. 3.1a, predicted modal amplitudes for the leading three harmonics are compared with hot wire measurements 
for the H18S12 case, where 18 micron roughness elements were applied at A, z = 12 mm. Effect of varying roughness 
height on the fundamental disturbance amplitudes is depicted in Fig. 3.1b. An analogous comparison for the 12 mm 
cases was presented in [10], wherein the secondary instability of the nonlinear stationary crossflow vortices was 
examined using temporal analysis. Because an envelope method was used to compute the integrated amplification 
of the secondary instabilities in [10], (unlike the present spatial analysis) the computed N-factors did not track the 
growth of a fixed disturbance entity. It should be noted that, whereas the unperturbed boundary layer flow in [10] 
had been predicted using the measured surface pressure distribution, the present set of results is based on a Cp 
distribution obtained by solving Euler equations for the conditions of interest. 

N-factor evolution for the unstable modes at various selected frequencies is shown in Figure 3.2 for each of the 
selected 12 mm cases. Corresponding growth rates for the 18-micron roughness (H18S12) case are also shown for 
illustration. Several secondary instability modes are found to exist in each of these cases, of which a few relevant 
modes were selected for plotting in Figure 3.2 and are represented by green, pink, and blue colors, respectively. The 
red curves denote lower frequency modes that originate significantly farther upstream as the traveling crossflow 
modes of the unperturbed boundary-layer flow but are substantially modified in the downstream region due to the 
finite amplitude stationary crossflow instability. 

The first set of (dominant) high-frequency secondary instability modes to become unstable is denoted by the green 
curves. It was determined from the relative contributions to the energy production terms that the green modes are 
associated with the wall-normal shear of the basic state (i.e. represent a y-mode in the notation of [10]) and attain a 
peak growth rate at around f=4.7 kHz for the 6- and 18-micron roughness cases and around f=4.9 kHz for the case of 
highest initial crossflow amplitude (i.e., 48 micron roughness case). This mode is denoted as y-mode. An additional 
set of modes (denoted via pink curves) becomes unstable either at about the same location as the y-mode, or 
somewhat farther downstream, depending on the roughness height used. This family of modes is associated with the 
span wise shear term and, hence, is denoted as z-mode 1 in this paper. It has a peak frequency of 1.6 kHz for lower 
initial amplitudes of the crossflow vortex and 2.9 kHz for the case with the highest initial crossflow amplitude. 
Subsequently, however, a second set of z-modes (denoted via blue curves) becomes unstable and achieves much 
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higher growth rates (and, hence, N-factors) in comparison with the first z-mode. The peak growth rate of this 
dominant set of z-modes (z-mode 2) is found to be near f=3.4 kHz for nearly all initial crossflow vortex amplitudes. 

It is seen that the maximum growth rate of the y-mode (denoted by the green curve) occurs between x/c = 0.3 and 
0.4, while the peak growth of the dominant z-mode (denoted by blue curve) occurs in bewteen x/c = 0.45 and 0.5. 
At the measured transition location of x/c=0.52 for the 6-micron roughness case, these y- and z-modes reach N- 
factors of approximately 10 and 7, respectively. As the roughness height (or, equivalently, the initial crossflow 
amplitude) is increased, each of these secondary instability modes become unstable at progressively upstream 
locations; however, the corresponding peak growth rates are found to decrease at the same time. The net result is 
that the transition N-factor attained by the y-mode 1 is increased to approximately 12, while that reached by the z- 
mode 2 is reduced to approximately 6 for the highest initial crossflow amplitude. N-factors for the intermediate 
roughness height (case H18S12) fall in between the corresponding values for the H6S12 and H48S12 cases, so that 
the overall spread in transition N-factors for either of the above two modes is small. The present N factors for the z- 
mode are somewhat lower in comparison with [10]; however, the precise cause behind this difference could not be 
ascertained. 

The above findings confirm the strong correlation between the measured transition location and the predicted 
secondary N-factor, regardless of whether the N-factor correlation is based on the amplification of the y-modes 
(which were found to be insignificant to transition in the related measurements by White and Saric [28]) or the z- 
modes of secondary instability (which were found to be an important catalyst for initiating the process of laminar 
breakdown in [28]). Detailed measurements in [28] showed the peak fluctuations associated with the z-mode to be 
slightly above 3kHz, which is consistent with the most amplified z-type secondary instability modes as predicted in 
the current analysis. The corresponding mode shapes for the dominant z-modes were also found to be rather close 
(compare, for instance, the mode shape of z-mode 2 at the bottom of Fig. 3.5 with Fig. 10 from ref. [28]). The 
reason(s) why the y-modes of secondary instability did not play an active role during transition in the experiment 
may be related to the receptivity characteristics pertaining to this mode; however, additional work is necessary to 
clarify those reasons. The recent work by Bonfigli and Klocker [26] appears to offer partial clues in this regard. 

The measurements in [28] indicated an additional spectral peak near 200 Hz that corresponds to traveling crossflow 
vortices modulated by the span wise variations associated with stationary crossflow mode [3, 26]. At the larger 
roughness height (i.e., larger initial stationary crossflow amplitude), the low frequency traveling disturbances 
attained significantly large amplitudes. Just prior to the laminar breakdown, the amplitudes of the low-frequency 
traveling modes were comparable to those of the z-modes of seconday instability, so it’s possible that both types of 
modes jointly contributed to the onset of breakdown at the above test condition. The behavior of the low-frequency 
traveling modes for A. z = 12 mm (cases H6S12 through H48S12) is investigated next. 

The low-frequency modes are unstable even in the absence of the stationary crossflow vortices. However, the 
stationary crossflow vortices, after they attain a finite amplitude, can significantly modify the behavior of these 
traveling modes. The red curves in Figure 3.2 show the N-factors and the growth rates of the traveling waves. Direct 
comparisons of the N-factors and growth rates of the modulated traveling waves for all three roughness heights are 
shown in Figure 3.3. These traveling waves start to become unstable approximately at X / C =0.04, and begin to 
show the effects of roughness height at approximately X I C =0.07, obviously because the presence of the growing 
stationary vortices is beginning to be felt at that stage. Up to x I C =0.2, the maximum N-factor for the traveling 
waves is lowest in the case of stationary crossflow vortices with the highest initial amplitude (blue set of curves in 
Fig. 3.3a). Based on the energy production budget, it was found that the traveling waves in the upstream region ( x/c 
< 0.15) are primarily associated with the wall-normal shear of the basic state (i.e., these modes correspond to y- 
modes). This is to be expected because they originate from a region in which the flow field is nearly uniform in the 
z-direction. 

The upstream peak in the growth rate curves (near x / C = 0.10— >0.12) corresponds to frequencies in the vicinity of 
150Hz. The initially dominant set of traveling waves have lower amplification rates in the region downstream of the 
of first lobe in the growth rate curves. Traveling modes near 460Hz attain their peak growth rates between x / C = 
0.25 and x / C = 0,36 depending on the initial amplitude of the stationary crossflow mode (i.e., the roughness height 
parameter). The second lobe in the growth rate curve from Fig. 3.3b is associated with the spanwise shear of the 
basic state (which is induced by the relatively strong stationary crossflow vortex). Thus, it is not surprising that the 
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peak of this second lobe shifts progressively upstream as the roughness height is increased from 6 micron (red 
curves in Fig. 3.3b) to 48 micron (blue curves in Fig. 3.3b). 

To illustrate the differences between the structures of traveling crossflow vortices in the upstream and downsteam 
regions, respectively, the mode shape evolution for a traveling mode at f=200Flz in the H06S12 case is examined 
next (Fig. 3.4). This particular traveling wave is seen to be a y-mode at small x/c and transitions to a z-mode beyond 
x/c= 0.25. The amplitudes of these waves are concentrated closer to the wall than the secondary instability waves, 
which ride on parts of the crosssflow vortices that are farther away from the wall as shown in Figure 3.5 for the 18 
micron case. The mode shape of the modulated traveling wave looks quite similar to the spatial distribution of low 
frequency oscillations observed in the experiment of White and Saric [28] (compare, for example, the mode shape 
shown at the bottom of Fig. 3.4 here with Fig. 7 of reference [28]). 

Because the computational resources required for secondary instability analysis can exceed the resources required 
for primary instability analysis by an order of magnitude or greater, it is worthwhile to examine if the onset of 
secondary instability may be correlated (in an empirical sense) with certain features of the primary flowfield. As a 
partial attempt towards that goal, Fig. 3.6 examines the spanwise averaged velocity profiles (obtained by using 
nonlinear PSE analysis) along the direction of the inviscid streamline at a few select streamwise locations for each 
of the three cases examined above. Fig. 3.6 suggests that the onset of the y-mode (high-frequency) secondary 
instability approximately coincides with the locations where the spanwise averaged velocity profile (aligned with the 
inviscid streamline at the boundary layer edge) first develops a visually distinct inflection point higher up in the 
boundary layer. For ease of comparison, those profiles (corresponding to x/c ~ 0.35, 0.30 and 0.25 for the 6, 18 and 
48 micron cases, respectively) have been highlighted in Fig. 3.6. Although not shown, the onset of this inflection 
point is also approximately correlated to the appearance of an additional, weak peak in the rms disturbance profiles 
associated with the stationary crossflow vortex. 

Figure 3.7 indicates the effect of 8 mm stationary crossflow waves on both traveling waves with the same spanwise 
wavelength and the high frequency secondary instability. Again, the initial amplitudes of the A, z = 8 mm mode have 
been chosen to achieve an approximate match between the overall predicted evolution of the r.m.s. primary 
amplitudes and the corresponding measured data for roughness heights of 6 micron and 48 micron, respectively. 
Carrillo et al. [27] report that the controlled excitation of stationary crossflow modes at a subcritical spacing delayed 
the onset of transition to x4;=0.8 1 at the lower roughness height (h=6 microns); however, the transition front moved 
back upstream to x/c=0.61 as the roughness height increased to h=48 microns (case F148S8). Results presented 
below attempt to provide some insight into these findings by examining the behavior of secondary instabilities in 
each of these cases. 

Figures 3.8a and 3.8b show the N-factors of the non-stationary modes with A, z = 8 mm for the above two cases. For 
h=6 microns (case H6S8), the traveling primary waves dominate, reaching an N-factor of approximately 6 at 
x/c=0.6, while the two secondary instability modes (a y-mode and a z-mode) barely show up during the region of 
significant primary amplitudes at A. z = 8 mm. At the larger initial crossflow amplitude corresponding to h=48 
microns, the growth of the traveling crossflow modes is significantly reduced (N=4 near x/c=0.6); but, the y- and z- 
modes of secondary instability reach significantly larger maximum N factors (N=6 and N=2.5, respectively) relative 
to the h=6 micron case. These computational findings suggest that there is little danger of premature transition due 
to the 8-mm control input when h=6microns (consistent with the large transition delay measured in the experiment), 
but there exists the possibility of an adverse impact on the transition location due to the relatively stronger secondary 
instability in the h=48 micron case (especially if the y-modes can get naturally excited at the smaller vortex 
wavelength). Whether or not this adverse effect may have been responsible for the upstream shift in the measured 
transition location (from x/c=0.81 for H6S8 case to x/c=0.61 for F148S8) cannot, however, be established on the 
basis of the present analysis alone. 


IV. Spatial Secondary Instabilities of Compressible Stationary Crossflow Vortices 

This section presents selected results for a Mach 2.4, 73-degree swept wing geometry that was modeled after the 
experimental configuration of Saric and Reed [14]. The nonlinear development of stationary crossflow vortices with 
various spanwise wavelengths and initial amplitudes were computed in Ref. [11] for a chord Reynolds number of 



Re c = 16 million. The dominant stationary crossflow mode in this case corresponds to a spanwise wavelength of 3 
mm. The first harmonic of this naturally dominant mode (1.5 mm spanwise wavelength) provides the simplest 
(although not optimal) form of control input to weaken the growth of the naturally dominant stationary modes. 
Determination of appropriate control input parameters (e.g., wavelength and amplitude of the control mode) is, of 
course, a major issue in designing an effective laminar flow control system. Specifically, the control input must be 
large enough to provide the desired control action (viz., sufficient modification of the basis state in order to induce 
the required stabilization of the “dangerous” modes), however, it cannot be excessively large as to precipitate 
premature transition (and/or reduce the extent of transition delay, as briefly alluded to in section III in the context of 
the experiments in [27]). As shown below, the secondary instability analysis (in conjunction with receptivity 
predictions [29]) may provide useful guidance to help select an appropriate range of control input amplitudes. 

The streamwise evolution of primary disturbance and mean flow distortion of the stationary mode with X z = 3 mm 
with an initial amplitude of Aj =10 7 is plotted in Fig. 4.1a, and the analogous evolution of the X z = 1.5 mm mode 
with several initial amplitudes is shown in Fig. 4.1(b). As discussed in Ref. [11], X z = 3 mm, the amplitude of the 
dominant mode rises quite rapidly along the chordwise direction, even when the amplitude is large enough to induce 
high-frequency secondary instability. The amplitude of the subcritical mode (Fig. 4.1b) reaches its peak 
significantly earlier than the A, z = 3 mm mode (Fig. 4.1a) and decreases rapidly thereafter. Thus, both the control 
action and the potential for premature transition due to such subcritical modes is limited to a finite spatial region in 
the vicinity of the modal peak. 

Amplification characteristics of secondary instability modes for the X z = 3 mm case are shown in Fig. 4.2. The N- 
factor curves for a range of selected frequencies are plotted in Fig. 4.2(a), whereas Fig. 4.2(b) shows the 
corresponding spatial growth rates. High-frequency secondary instability modes are seen to exist for frequencies up 
to 2000 kHz. The maximum N-factor for the secondary modes reaches 10 even before the peak primary amplitude 
has been reached, indicating that (if, indeed, the primary instability spectrum were to mirror the inflow behavior in 
Fig. 4.2(a) and the disturbance environment can excite the relevant secondary instability modes) the onset of 
transition should occur near x/c = 0.6. The growth characteristics of the secondary modes also support the 
hypothesis [11] that, in this particular case, movements in transition onset might approximately correlate with the 
corresponding shift in the region of rapid rise in the primary amplitude. 

A more detailed examination of the X z = 3 mm case reveals that there exist at least four secondary instability modes, 
all of which are y-modes and two are shown in Fig. 4.2 (b). A y-mode is the first to become unstable at 
approximately x/c = 0.5, this mode reaches an N-factor of approximately 8.5 for a frequency of approximately 1.05 
MHz at x/c = 0.60. However, transition is likely to be caused by a second y-mode with a frequency near 1.3 MHz, 
which becomes unstable at x/c ~ 0.55, overtakes the y-mode 1 in terms of N-factor near x/c = 0.57 and reaches a 
maximum N-factor of approximately 10 at x/c = 0.60. The other two y-modes (not shown) reach much lower N- 
factors. 

The next set of results pertains to the variation in secondary instability N-factors as the initial amplitude of the 
control mode (A, z =1.5 mm) is varied (assuming this mode alone to be dominant over the relevant spatial region). 
The effect of modulation of primary instability by stationary crossflow vortices of initial amplitudes of 0.001, 0.002 
and 0.005 is shown in Figure 4.3, and secondary instability results for initial amplitudes of A, = 0.001 and 0.002 
(which yields a modest stabilization of the X z = 3 mm stationary mode as described in [11]) are shown in Figs. 
4.4(a)-(b) and 4.4(c)-(d), respectively. 

Figure 4.3 shows that the growth rates of the modulated primary instability immersed in higher amplitude crossflow 
vortices have lower growth rates up to x/c ~ 0. 14, indicating the weak stabilizing influence of moderate amplitude 
stationary crossflow mode on lower amplitude traveling crossflow modes. Farther downstream, however, this trend 
is reversed. Higher amplitude stationary crossflow vortices now appear to enhance the linear growth of the traveling 
modes that are dominant in the downstream region, as reflected by the progressively upstream shift in the (local) 
peak in the growth rate curves as the roughness height is increased (similar to the low-speed case examined in 
section III). However, the linear growth mechanism takes on different characteristics as the growing traveling modes 
can now be either z-modes and/or y-modes. 
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For X 7 = 1.5 mm, up to six different high-frequency secondary instability modes were found, with three predominant 
ones (Figure 4.4). Two of those three modes are y-modes and the remaining one is a z-mode. They are shown via 
different colors in Figure 4.4. For the A, = 0.001 case, the three modes have peak frequencies of 0.88, 1.72 and 0.72 
MHz, respectively, and for the A, = 0.002 case, those frequencies increase to 1.16, 2.0 and 1.56 MHz, respectively. 
For an initial amplitude of 0.001, the maximum N-factor approaches 9 for one of the y-modes. For an initial 
amplitude of 0.002, the N-factor reaches approximately 13. All of these findings indicate a strong possibility of 
premature transition at this level of control input. The finding that strong secondary instabilities can occur at both 
A, z = 3 mm and K, = 1.5 mm (and, therefore, at the intermediate wavelengths as well) is suggestive of the potentially 
delicate nature of roughness based transition control at higher chord Reynolds numbers. 

The progressive distortion of the mode shape corresponding to the (primary) traveling crossflow modes over this 
supersonic configuration is indicated in Fig. 4.5. The distorted mode shape of the traveling crossflow mode at 
x/c=0.25 may be compared with the velocity contours of the base flow (top of Fig. 4.6) and the local mode shapes of 
the dominant secondary instability modes (upper and lower middle, and bottom of Fig. 4. 6). The associated energy 
production budget reveals that both the wall-normal and spanwise shears of the basic state make comparable 
contributions to the growth rate of this traveling crossflow mode at x4.~0.25, 

A necessary condition for a 2D compressible shear flow to develop inviscid type of instability is the existence of a 
generalized inflection point. To assess whether a significant change in the inflection point characteristics of the 
spanwise averaged velocity profile may serve as a crude indicator for the onset of secondary instability, the density 
weighed shear profiles for the velocity component along the inviscid streamline direction are now examined. The 
<p5U s /3y> profile in the supersonic case exhibits a local maximum corresponding to a generalized inflection, which 
is present at all locations regardless of the primary disturbance amplitude. However, after the primary amplitude 
becomes sufficiently large, the profile develops an additional unstable inflection point. For the 3 mm cases the 
appearance of this additional generalized inflection point is close to the location where the dominant secondary 
instability mode first becomes unstable. For the 1.5 mm cases, the generalized inflections occur at x/c= 0.19 and 
0.17, respectively, for the two initial amplitudes considered herein (Fig. 4.7). In these cases, the onset of inflection 
is somewhat ahead of the onset of the y-type secondary instability. 


V. Summary and Concluding Remarks 

This paper presented planar (2D) eigenvalue analysis methodology for spatial instabilities of compressible shear 
flows with two inhomogeneous directions, such as crossflow dominated swept-airfoil boundary layers. Direct 
spatial 2-D eigenvalue computations, rather than approximate calculations based on a combination of temporal 
analysis and the well-known Gaster’s relation, were performed to characterize the growth characteristics of 
secondary instabilities of stationary crossflow vortices. The equations governing these instabilities were formulated 
in a non-orthogonal coordinate system, so that the proper periodic boundary conditions may be imposed in the 
direction parallel to the leading edge of an infinite-span swept airfoil and the exponential growth in the direction 
perpendicular to the leading edge can be computed without any ambiguity. 

The choice of numerical case studies was motivated by transition prediction and control for swept wing boundary 
layers. The selected swept-airfoil configurations included (i) the low-speed, 45-degree sweep, NLF-0415(b) 
configuration with a chord Reynolds of Re c = 2.4 million, as tested previously in the ASU Unsteady Wind Tunnel, 
and (ii) a Mach 2.4, 73-degree sweep configuration with Re c = 16 million. For each of these configurations, we 
examined both the modification of the traveling crossflow modes and the destabilization of high-frequency 
secondary instabilities in the presence of finite amplitude stationary crossflow vortices of a specified spanwise 
wavelength and varying initial amplitudes. 

Regardless of the speed and/or the chord Reynolds number, the onset of secondary instability was found to move 
forward as the initial amplitude of the stationary crossflow vortex was increased. As the stationary vortex increased 
in amplitude along the downstream direction, the base flow along the direction of the inviscid streamline was shown 
to develop a inflectional (or, in the compressible case, a generalized inflectional) character higher up in the boundary 
layer and the locations for the onset of the high-frequency secondary instability seemed to fall at, or somewhat ahead 
of , the emergence of this inflection point. 
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The stationary crossflow vortices also modulate any small-amplitude traveling crossflow vortices, resulting in a 
weak stabilization of the initially dominant traveling modes. The initial stabilizing effect is more pronounced for 
higher initial amplitudes of the stationary crossflow vortex. Farther downstream, however, the above trend is 
reversed such that the most unstable traveling modes in that region are actually destabilized as a result of the 
stationary mode. 

Overall, the results reconfirm that secondary instability analysis can explain a number of observed features during 
the previous experiments in the ASU Unsteady Wind Tunnel, both with and without the use of surface roughness at 
a subdominant wavelength as a means of transition control. Application of the 2D eigenvalue analysis to the higher 
Reynolds number, supersonic configuration revealed that a broad spectrum of stationary crossflow modes (including 
modes that would be suitable for roughness based transition control) can sustain sufficiently strong secondary 
instabilities, so as to potentially induce transition over this configuration. 

Acknowledgments 

The authors are grateful to Dr C.-L. Chang for his assistance and for useful technical discussions. The present work 
was supported by the Supersonic Cruise Efficiency (Airframe) discipline of the Supersonics Project under NASA's 
Fundamental Aeronautics Program, with Ms. Linda Bangert as the Associate Principal Investigator for this 
discipline. 

References 

1. Pfenninger, W. and Bacon, J. W., “Amplified Laminar Boundary Layer Oscillations and Transition at the Front 

Attachment Line of a 45 ° Flat-nosed wing with and without boundary Layer Suction”, Viscous Drag 
Reduction, ed. C. S. Wells, Plenum, 1969. 

2. Mack, L. M. “Boundary Layer Stability Theory”, Report 900-277, Rev. A, Jet Propulsion Laboratory, Pasadena, 
CA, November 1969. 

3. Malik, M. R., Li, F. and Chang C.-L., "Crossflow Disturbances in Three-Dimensional Boundary Layers: Non- 
linear Development, Wave Interaction and Secondary Instability." J. Fluid Mech. Vol. 268, 1994, pp. 1-36. 

4. Bushnell, D.M. and Malik, M.R., “Application of Stability Theory to Laminar Flow Control - Progress and 
Requirements,” Stability of Time Dependent and Spatially Varying Flows,” D.L. Dwoyer & M.Y. Hussaini, 
eds.. Springer Verlag, pp. 1-17, 1987. 

5. Malik, M. R., "Hypersonic Flight Transition Data Analysis Using Parabolized Stability Equations with 
Chemistry Effects ” J. of Spacecrfat and Rockets 40,2003, pp. 332-344. 

6. Reibert, M. ,Saric , S. Carrillo, R.B., Jr and Chapman, K. L. “Experiments in Nonlinear Saturation of Stationary 
Crossflow Vortices in a Swept- wing Boundary Layer,” AIAA Paper 96-0184, 1996. 

7. Kohama, Y. Ondodera, T. and Egami Y. “Design and Control of Crossflow Instability Field,” In Proc. IUTAM 
Sym. On Nonlinear Instability and Transition in Three Dimensional Boundary Layers,” Manchester, U.K (ed. 
P.W. Duck and P. Hall). Pp- 147-156. Kluwer, 1996 

8. Saric, W.S., Reed H.L., and White E.B. “Stability and Transition of Three-Dimensional Boundary Layers,” 
Ann. Rev. Fluid Mech., vol. 35, pp. 413-440, 2003. 

9. Saric, W.S., “Gortler Vortices,” Ann. Rev. Fluid Mech., vol. 26, pp. 379-409, 1994. 

10. Malik, M. R., Li, F. Choudhari, M. M. and Chang C.-L., "Secondary Instability of 
Crossflow Vortices and Swept-wing Boundary-layer Transition." J. Fluid Mech. Vol. 399, 1999, pp. 85- 1 15. 

11. Choudhari, M., Chang C.-L, Streett C.L., and Balakumar, P., “Integrated Transition Prediction: A Case Study in 
Supersonic Laminar Flow Control,” AIAA Paper 2003-0973 , 2003. 

12. Choudhari, M, Chang, C.-L., and Jiang, L., "Towards Transition Modeling for Supersonic Laminar Flow 
Control,” Philosophical Transactions of Royal Society of London (Physical and Mathematical Sciences), vol. 
363, no. 1830, pp. 1079-1096, May 2005. 

13. El-Hady, N.M. “Spatial three-dimensional secondary instability of compressible boundary-layer flows,” AIAA 
Journal, vol. 29, no. 5, pp. 688-696, 1991. 

14. Saric, W. S. and Reed, H. L., “Supersonic laminar flow control on swept wings using distributed roughness.” 
AIAA paper 2002-147, 2002. 

15. Hall, P., and Horseman, N., "The Linear Inviscid Secondary Instability of Longitudinal Vortex Structures in 
Boundary Layers, ” J. Fluid Mech., vol. 232, pp. 357-391, 1991. 

16. Chang, C.-L., Malik, M. R. and Vinh, H. “Linear and Nonlinear Instability of Compressible Swept-Wing 
Boundary Layers.” AIAA Paper 95-2278, 1995. 


11 



17. Li, F. and Malik, M.R., “Fundamental and Subharmonic Secondary Instability of Gortler Vortices,” , J. Fluid 
Mech. vol 297, pp.77-100, 1995 

18. Lin, R.-S., Malik, M. R., "On the Stability of Attachment-line boundary Layers 

Part 1. The Incompressibbe Swept Fliemenz Flow." J. Fluid Mech. Vol. 311, 1996, pp. 239-255. 

19. Lin, R.-S., Malik, M. R., “Stability and Transition in Compressible Attachment-Line Boundary-Layer Flow”, 
SAE Technical Paper No. 952041, Aerotech ’95, Los Angles, CA, September 18-21, 1995. 

20. Theofilis, V., “Advances in global linear instability analysis of nonparallel and three-dimensional flows,” 
Progress in Aerospace Sciences, vol. 39, no. 4, pp. 249-315, May 2003. 

21. Janke, E. and Balakumar, P., “On the secondary Instability of Three-Dimensional Boundary Layers,” Theoret. 
Comput. Fluid Dynamics, vol. 14, ppl67-194, 2000. 

22. L. Jiang, M. Choudhari, C.-L. Chang, and C. Liu, “Direct Simulation of Instability-Wave Generation and 
Propagation in Supersonic Boundary Layers,” in Proceedings of International Conference on Computational 
Science and Its Applications - ICCSA 2003, Montreal, Canada, May 18-21, 2003, Vol. II (Lecture Notes in 
Computer Science 2668), V. Kumar, M.L. Gavrilova, C.J. K. Tan, and P. LEcuyer (Eds.), Springer Verlag, 
Berlin-Heidelberg, July 2003. 

23. L. Jiang, M. Choudhari, C.-L. Chang, and C. Liu, “Direct Numerical Simulation of Crossflow Disturbances in 
Supersonic Boundary Layer,” AIAA Paper 2004-0589, 2004. 

24. Chang, C.-L., “The Langley Stability and Transition Analysis Codes (LASTRAC: LST, Lineara dn Nonlinear 
PSE for 2D, Aximsymetric, and Infinite Swept Winga Boundary Layers”, AIAA Papert 2003-0974, 2003. 

25. Reibert, M. Saric, W. S., Carrilo, R. B., Jr & Chapman, K. L. 1996 Experiments in Nonlinear Saturation of 
Stationary Crossflow Vortices a a Swept Wing Boundary Layer. AIAA paper 96-0184. 

26. Bonfigli, G, and Kloker, M., “Secondary instability of crossflow vortices: validation of the stability theory by 
direct numerical simulation,” J. Fluid Mech., Vol. 583, pp. 229-272, 2007. 

27. Carrillo, R. B., Jr, Reibert, M. S. & Saric, W. S. 1996 “Distributed Roughness Effects on Stability and 
Transition in Swept Wing Boundary Layers,” Final report for NASA LaRC co-op. agreement NCC- 1-194. 

28. White E. B., and Saric W. S., “Secondary Instabilty of Crossflow Vortices.” J. Fluid Mech. Vol. 525, pp. 275- 
308, 2005. 

29. Chang, C.-L., and Choudhari, M, “Boundary-Layer Receptivity and Integrated Transition Prediction,” AIAA 
Paper 2005-0526, 2005. 



Figure 2.1. Schematic of a swept wing of infinite span, with three surface coordinate systems, namely, 
(x, Z.) , (X , Z) and (£,£), the last of which corresponds to the noil-orthogonal, vortex-aligned coordinate 
system. 
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Figure 3.1 (a) Comparison of experimental and 
computational modal amplitudes for 12 mm crossflow 
vortices in case H18S12. 



Figure 3.1 (b) Comparison of measured and predicted 
fundamental amplitudes for cases H06S12, H18S12, 
and H48S12 (curves: computation, symbols: 

experiment). 




Figure 3.2 Top left, top right and bottom left: N-factor curves for selected modes for 6-, 18- and 48-micron 
roughness cases, respectively. Bottom right: spatial growths rates for the 18-micron case. Measured transition 
locations are indicated by a dashed red line. Color of N-factor curves denotes the mode type, whereas different 
curves for a given mode correspond to disturbances at different frequencies. The most amplified frequency for 
each selected mode is highlighted on the plot. 
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x/c 


Figure 3.3 (a) N -factors for modulated traveling 
crossflow vortices for different roughess heights. Color 
of curve indcates roughness height, whereas different 
curves of the same color are for different frequencies. 



Figure 3.4 Eigenfunction evolution for modulated 
traveling mode at f=200 Hz. The abscissa and ordinate 
correspond to suitably normalized spanwise coordinate 
£ and wall-normal coordinate y, respectively. 



Figure 3.3 (b) Growth rates of modulated traveling 
crossflow vortices corresponding to cases in Fig. 
3.3(a). Peak growth-rate frequency in each case is 
indicated. 



Figure 3.5 Magnitude of secondary instability 
eigenfunctions for velocity component along the vortex 
axis for 45-deg swept NLF-0415(b) wing. Top: y-mode, 
middle: z-mode 1 and bottom: z-mode 2 (Abscissa and 
ordinate similar to Fig. 3.4). 



Figure 3.6. Mean velocity profiles aligned with the local inviscid streamline. Profiles begin at x/c=0.05 and end 
at x/c=0.6, with an increment of 0.05 across successive profiles. Highlighted profile locations indicate where the 
profile develops a distinctly inflectional character. 


14 





Figure 3. 7. Comparison of experimental and computational r.m.s. amplitudes of stationary disturbance for a 
roughness array spacing of 3^=8 mm (cases H06S08 and H48S08 for the NLF0415(b) configuration). 


H06S08 


Modulated primary mode 

z-mode 

y-mode 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

x/c 



H48S08 

Modulated primary mode 

y-mode 1 4.2 kHz 

z-mode 

y-mode 2 , 1 .7 kHz 


2.1 kHz 


0.21 kHz 


(a) 6-micron roughness (b) 48-micron roughness. 

Figure 3.8. N -factor curves for 3^=8 mm ( dashed line indicates measured transition location; line color denotes 
selected mode type, whereas different curves of same color indicate different frequencies for a given mode; 
highlighted frequencies correspond, approximately, to most amplified frequency for each selected mode). 


15 




(b) X, = 1. 5mm for initial fundamental amplitudes 


(a) X z = 3mm 


of 0. 001 and 0. 002 

Figure 4.1. Primary disturbance evolution along Mach 2.4 swept-wing boundary layer. Amplitudes of several 
disturbance harmonics, including the fundamental mode ( mode 1) and mean-flow-correction (mode 0), are 
shown in each figure. 




Figure 4.2(a) N-factor curves. Figure 4.2(b) Spatial growth rate curves. 

Figure 4.2. Secondary instability of stationary crossflow mode at X — 3mm (At = le-7 in Fig. 4.1(a)). Line 
color denotes selected mode type, whereas different curves of same color indicate different frequencies for that 
mode; highlighted frequencies correspond, approximately, to most amplified frequency for each selected mode. 
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Figure 4.3. Effect of initial amplitude of the stationary crossflow vortex on growth rates of modulated 
traveling crossflow instability at A. = 1.5 mm (Line color denotes initial amplitude. Red curve: A,- = 0.001, 
Green: Ai = 0.002, and Blue: A; = 0.005. Different curves at each initial amplitude represent different 
frequencies, with peak growth rate frequencies indicated on the plot). 
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Figure 4.4(b) Spatial growth rate curve, A ; = 0.001 
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Figure 4.4(c) N -factors curves, A, = 0.002 
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Figure 4.4(d) Spatial growth rate curves, A,- = 0.002. 


Figure 4.4. Secondary instability of stationary crossflow mode at A_. = 1.5 mm. (Different curves of the same 
color are for different frequencies, and peak growth rate frequencies for each mode are indicated in parts b and 
d). 
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Figure 4.5. Eigenfunction contours for modulated Figure 4.6. Basic state and secondary mode shapes at 
traveling crossflow waves at 120 kHz, (A, = l .5 mm ; x/c = 0.25 ( A_ — 1 .5 mm. ). Abscissa and ordinate 


x/c ~ 0.08, 0.11, 0.16, and 0.25 from top to bottom, 
respectively). Abscissa and ordinate similar to Fig. 
3.4; Contour variable corresponds to magnitude of 
eigenfunction for the velocity component along the 
direction of stationary vortex ( which is approximately 
aligned with the local inviscid streamline ). 


similar to Fig. 3.4. Top: Contours of velocity 
component along the axis of the stationary crossflow 
vortex; Upper and lower middle, and bottom: 
Secondary instability eigenfunction corresponding to 
velocity component along the vortex axis for two y- 
rnodes and one z-mode, respectively, at f =1700 kHz. 
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Figure 4.7 Spanwise averaged p3U s /3y profiles aligned with the direction of the inviscid streamline for the 
supersonic swept-airfoil configuration (X,=3 and 1.5 mm), with highlighted locations indicating where the profile 
first develops an additional generalized inflection due to the finite amplitude of stationary crossflow vortices. 
Profiles start at x/c =0.1 and proceed in intervals of Ax/c = 0.05, with the exception of x/c=0.54 (left plot), 
x/c=0.19(middle plot) and x/c=0.17 (right plot). 
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